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Abstract: Using the Lewis-Riesenfeld method of invariants we construct explicit an¬ 
alytical solutions for the massless Dirac equation in 2-1-1 dimensions describing quasi¬ 
particles in graphene. The Hamiltonian of the system considered contains some explicit 
time-dependence in addition to one resulting from being minimally coupled to a time- 
dependent vector potential. The eigenvalue equations for the two spinor components 
of the Lewis-Riesenfeld invariant are found to decouple into a pair of supersymmetric 
invariants in a similar fashion as the known decoupling for the time-independent Dirac 
Hamiltonians. 


The two dimensional massless Dirac equation has recently attracted a lot of renewed 
attention because it describes quasi-particles in graphene |l], which is well known 

to possess a large amount of remarkable properties. Especially the Dirac equation in the 
presence of a magnetic field is of great interest as, unlike electrostatic potentials, such 
a conhguration allows in principle to conhne the Dirac fermions i i I- Many exact 
solutions have been provided for a variety of time-independent Hamiltonians and mag¬ 
netic field configurations ^ including some for complex magnetic helds leading to 
pseudo/quasi-Hermitian interactions ©0- While some solutions for the time-dependent 
Dirac equation in 1-|-1 dimensions have been constructed |12, |^, 14], little is known about 
the time-dependent setting with a magnetic field in 2-|-l dimensions and no exact solutions 
have been reported. The aim of this manuscript is to commence filling that apparent gap. 
We shall demonstrate that the Lewis-Riesenfeld method of invariants [IT] is a technique 
which can be employed successfully to solve this problem. 

We consider here the time-dependent massless Dirac equation in two spacial dimensions 
in the form 

H'i' = idt'^, with H{x,y,t) = a ■ p, (1) 


and the two component wave function T' = The effective Hamiltonian includes 

comprised of the standard Pauli matrices and p being the two-dimensional 
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momentum vector minimally coupled to a vector potential A = (Ax, Ay) in the standard 
fashion p = (a(t)px + Ax,b(t)py + Ay). Besides the time-dependence entering through 
the vector potential, we allow here also for explicit time-dependent factors in front of 
the momenta, a{t) and h{t). One possibility to think of these factors is that they result 
from a time-dependent background, as discussed in |]^] or alternatively as time-dependent 
velocities. The velocity of light, the reduced Planck constant and the charge are all set to 
one in our discussion, i.e., c = h = e = 1. 

We make now some specially simply choices for the vector potential by taking Ax = 0 
and Ay(x,t) = g{t)x. In addition, we assume that the wave function separates with a free 
plane wave moving in the y-direction 

Af(x,y,t)=e^^’^y-‘^^^^ix,t), (2) 

with wave number k and frequency cu, such that we are left with the task to solve 

= idt^, with H(x,t) = a{t)axPx + [kb{t) g(t)x] ay — wl, (3) 

for the two component wave function $ = (4>i,(l)2)- The first observation we make here is 
that the form of the explicit time-dependence does not allow for the standard decoupling of 
the systems into a pair of Hamiltonians related to each other by intertwining operators as 
common in the time-independent Dirac equation in analogy to standard supersymmetric 
quantum mechanics ||T^, |l^ . 

We will attempt to solve equation (^) by using the Lewis-Riesenfeld method orig¬ 
inally designed to solve the time-dependent Schrodinger equation. The first step in this 
approach consists of solving the evolution equation 

^=8i/(t) + j|/(t).H(«)l=0, (4) 

for the Hermitian time-dependent invariant I(t). As usual in this context we take the 
invariant to be of the same order and form in the canonical variables as the Hamiltonian 

I{t) = a{t)pxP(t)x-\-'y{t), (5) 

where a(t), /3(t) and 7 (t) are now unknown time-dependent matrices. Substituting our 
Ansatz into the evolution equation yields six constraining equations for the three 
coefficient matrices 


[a,ax]=0, [P,ay] = 0, g[a,ay]a[^,ax] = 0, ( 6 ) 

—ia = kb [a, Uy] -)- a [ 7 , Ux], (7) 

-i'l3 = g['y,ay\, ( 8 ) 

—^7 = kb [ 7 , ay] + iaax/3 — igaay. (9) 


Expanding the matrices in the su(2)-basis, a(t) = ai(t)I+a 2 {t)ax + a 3 (t)ay + ai(t)az with 
ai(t) G M for f = 1 , 2 ,3,4 and similarly for (3{t), 7 (f), these equations are straightforward 


- 2 - 




Time-dependent massless Dirac fermions in graphene 


to solve. Starting with the first two equations immediately imply that as = 04 = 0 
and /32 — P 4 — 0- The last equation in (^) then yields /3^ = a 2 g{t)/a{t). Proceeding in 
this way for we find the following form for the time-dependent invariant 

I(t) = {aipx 7i) I a2PxCrx + [fd^x + 73(t)] ay, (10) 

with constants 01 ,^ 1 , 02 , The time-dependence of I{t) is entirely contained in the 
function 73 (t) which is constrained by 

73 W = aig{t), and = ka 2 ^ = kjdj^. (11) 

a{t) g{t) 

In addition we found that a{t) = pig{t) with p, = a 2 /(d^ being constant has to be satisfied. 
The equations ( p(T|) are most conveniently solved in terms of b{t) 

!■{»)*) ^ . 9{t) = 

The next step in the Lewis Riesenfeld approach consists of solving the eigenvalue 
equation for the time-dependent invariant, i.e., we need to solve 


73W = J 



= Xx{t), 


(13) 


for the time-dependent eigenfunction x(t) = (x_|_(t), x_(t)) and time-independent eigen¬ 
values A. For this purpose we note at first that we can write ([I^ as 


0 

T+ 




— (A-ai - 7 i) 
02 



(14) 


Thus we notice that unlike as the time-dependent Hamiltonian the invariant equation can 
be decoupled easily and acquires the form of a supersymmetric pair. Acting again with the 
off-diagonal invariant operator on (U) we obtain the two decoupled equations 

I±x± = [pI + W^± it'] x± = - «1 - 7 i)'x±, (15) 

"2 


for the two operators I± := L^L±, where W{x,t) = [ 73 (t) -|- Id^x]/a 2 is the analogue to 
the superpotential in standard supersymmetric time-independent quantum mechanics. We 
observe that the potential is still time-dependent, but now a simple re-definition of our 
variables will move this dependence entirely into x±- Defining the new time-dependent 
variable f,{t) = [ 73 (t) -|- (d^x]/a 2 converts (|l^ into two eigenvalue equations for the time- 
independent quantum harmonic oscillator 


1 ^ 

'2df 


(-2 
2 ^ 


± 


X± = ^ X±: 


(16) 


with 2A± = (A — ai — 7 i )^//33 =F P- Demanding x±(.^) to be a square integrable function 
L^(M, d^) the solution is of course 

X±,n(0 = O ^ = P , (17) 
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with Hn denoting the n-th Hermite polynomial. This means the two eigenvalues for the 
time-dependent spinor components of the invariant in (13) quantize to 

= ai + 7 i + |33S^y^J,{2n + l±l), (18) 

with s = ±1 being two possible signs of the square root. As expected these eigenvalues are 
indeed time-independent. Due to the supersymmetric structure we have the standard shift 
in the eigenvalues, that is 

As argued by Lewis and Riesenfeld [R| the eigenfunction of the Hamiltonian and the 
invariant just differ by a phase 

= (19) 


where the real function 5{t) in (19) must obey 

d6{t) 


dt 


= {Xn\idt- H{t)\Xn) 


( 20 ) 


The right hand side of (2C) can be computed directly with our known eigenfunctions. We 
obtain 


{Xn\idt - H{t)\Xn) = {x+,n\idt + uj\x+,n) + {x+,n\ - aPx + i{kbgx) \x-,n) (21) 

-h {x-,n\ - apx - i{kb + gx) \x+,n) + (x-,n| + OJ\x-,n) 

= 2UJ, ( 22 ) 

where we used the orthogonality of the eigenfunctions, dt |x±,n) = 0) {x+,n \ Px \ x-,n) 

- {x-,n\px\x+,n) and (x+,n| X\x-,n) = {X-,n\^\x+,n)- The phase therefore simply be¬ 
comes 

6{t) = 2ut. (23) 

Thus assembling the results from equations dD, (0), (H) and provides an exact 
solution to the Dirac equation (||). 

We have demonstrated that the Lewis-Riesenfeld method can be applied to construct 
solutions to the 2-|-l dimensional time-dependent Dirac equation. The time-dependence 
resulted from a background and a magnetic field. Exact solutions for this type of scenario 
have not been known previously. Clearly there are plenty of open problems and challenges 
left. For instance, just as in the time-independent scenario one would like to know exact 
solutions for more complicated vector field configurations, different background scenarios 
and possibly different assumptions about the motion in the y-direction. These tasks are 
left for future work, where this note can be taken as encouragement as it demonstrates the 
successful applications of a method to tackle these kind of problems. 
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